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Abstract 

Initial-boundary value problem for linearized equations of motion of 
viscous barotropic fluid in a bounded domain is considered. Existence, 
uniqueness and estimates of weak solutions to this problem are derived. 
Convergence of the solutions towards the incompressible limit when com- 
pressibility tends to zero is studied. 

1 Introduction 

In many cases mathematical treatment of liquids is done in the framework of 
incompressible fluid. However, from the physical point of view, all the liquids 
existing in nature are low compressible. Therefore it is reasonable to study 
asymptotic properties of solutions to equations of low compressible fluid, in 
particular, convergence to the corresponding incompressible limit. 

Low Mach number limit, which can be considered as a particular case of 
low compressibility limit, was studied by E. Feireisl, P.-L. Lions and others 
[TJ [31 [3111]. In particular, it was proved (see that there exists a sequence 

of weak solutions to equations of compressible fluid motion which converges to 
a solution of the corresponding equations of incompressible fluid. More pre- 
cisely, weak convergence of the velocity was established. But from the physical 
consideration, one could desire strong convergence of the solutions to yield bet- 
ter approximation of compressible fluid by incompressible one. Therefore it is 
interesting to study sufficient conditions for the convergence to be strong. 

Strong convergence of the velocity was established in [5j III, §8] for the solu- 
tions of "compressible" system arising in the method of artificial compressibility. 
It was also proved that the gradient of the pressure converges weakly, but strong 
convergence of the pressure was not examined. 

In some situations the convergence to the incompressible limit cannot be 
strong. In |^ there was derived a necessary condition of strong convergence 
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of classical solutions to equations of compressible fluid when compressibility 
tends to zero. The condition represents a restriction on the initial condition for 
the equations of incompressible fluid and this restriction cannot be satisfied in 
general case. 

In this paper we study weak solutions to linearized equations of compress- 
ible fiuid. The reason for this is twofold. First, by the present moment only 
existence of weak solutions (on arbitrary time interval) to the full equations of 
compressible fluid is established [2] , and there is well-known regularity problem 
for these solutions. Second, the proper treatment of nonlinear terms not only 
would be hard technically but might perhaps obscure the effect of low compress- 
ibility on the solutions, which is the main subject of this paper. (Also note that 
the linearized equations of compressible fluid are of interest on they own. For 
instance, spectral properties of the operator corresponding to linearized steady 
equations of compressible fluid were examined in [7].) 

Linearized equations of compressible fluid motion were studied in [51 [S]- 
Estimate of strong generalized solution to initial-boundary value problem for 
these equations in a bounded three-dimensional domain was derived in [S] , and 
existence of strong generalized solution to Cauchy problem in the whole space 
for them was established in [5]. It appears that existence of weak solutions to 
these equations has not been addressed. 

In this paper we derive existence, uniqueness and estimates of weak solutions 
to the initial-boundary value problem for linearized equations of compressible 
fluid. We examine convergence of the solutions to the incompressible limit when 
compressibility tends to zero. Briefly, we prove that 

• in general case the velocity field converges weakly; 

• if initial condition for the velocity is solenoidal then the velocity converges 
strongly and the pressure converges weakly; 

• if, in addition, the initial condition for the pressure is compatible with the 
initial value of the pressure in the incompressible fiuid then the conver- 
gence of the pressure is strong. 

2 Notation and Preliminaries 

2.1 Common Functional Spaces 

Let E C M" be a bounded domain n e N. We will use the following standard 
spaces: 

• LP{E) is the Lebesgue space of real-valued functions on E, 1 < p < oo; 

• LP{E) = {u\ue LP{E), Jj^udx = 0}; 

• H^i^E) = W^'P{E) is the Sobolev space of real- valued functions whose 
weak derivatives up to order s e N belong to LP{E); 
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• D'{E) is the space of distributions on E; 'D{E) is the space of test functions 
on E; 

• C^{E) is the space of smooth real- valued functions on E with compact 
support; 

• H^{E) denotes the closure of C(f (S) in H\E)-noYm; 

• H^^{E) denotes the dual space of Hq{E); 

(For E — {a,b) ClR we wiU omit undue brackets, i.e. L^{a, b) — L^{{a, 6)).) 

For R'^-valued functions (k e N) we will use Cartesian products of these 
spaces, e.g. L'^{E)^ . Since H~^{D)'' is linearly and continuously isomorphic to 
the dual space of HI{DY, let H-^{Df denote the latter. 

Let C K"^ be a bounded domain with a partially-smooth boundary dD^ 
deN, d> 2. Let Qt = Dx (0, T), where T > 0. 

Vector-valued functions will be denoted by bold letters. For such functions 
we will use the following spaces: 

• j{D) = {u| u e C^{Df, divu ^ 0}; 

• J{D) is the closure in L'^{D)-nQTUY of j{D) (cf. 
(This space is often also denoted as H{D).) 

• V{D) is the closure in Hl{DY-noviTY of j{D). 
2.2 Scalar Products and Duality 

It is well-known that the space Hq(D)'^ (fc e N) is a Hilbert space with respect 
to the dot product 

((u,v))= / (V u) : (V (8) v) da; = / {d^Uk){diVk) dx 

J D J D 

where {u,v} C H}^{D)^ are vector- valued functions (here and further Einstein 
summation convention is used) . 

Let (•,•), depending on the context, denote the standard dot product in M*' 
or the dot product in L^{D)^, i.e. (u, v) = Jj^UiVidx. 

Let II • ||x denote the norm of a Banach space X (with dual space X*) and let 
(•, •)^ denote the duality brackets for the pair {X*,X). Let (a;„) C X , x E X . 
It is convenient to use the following notation: 

• Xn ^ X means that the sequence (a;„) converges to x weakly; 

• Xn —r X means that the sequence (x„) converges to x *-weakly. 

We will also use the following short-hand notation for norms (see [S]): 

I • I = II • llL2(£))fc 

II ■ II = II • ll-ffoM-D)* 

where the value of fc € N is defined by the context. 
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2.3 Differential Operators 



It is known that the operators V, A and div , defined on smooth functions with 
compact support (i.e. functions from C^{D)'^, k S N), can be extended by 
continuity to bounded Unear operators V: L^{D) H^^{D)'^, A: Hq{D)'^ 
H-^{DY^ and div : L'^{DY H~^{D), defined by 



2.4 Spaces of Banach Space Valued Functions 

Consider an arbitrary closed interval [0,T] where T > 0. Let X be a Ba- 
nach space and s e N. Let LP{0,T]X) and W'P(0,r;X) denote accordingly 
Lebesgue-Bochner and Sobolev-Bochner spaces of X- valued functions of real 
variable t G [0, T] (see, e.g., HJ), 1 < p < oo. Let ft denote the weak derivative 
of a function / e W^'P{0, T; X) with respect to t. 

Now we recall the following well-known property of Sobolev space of Banach 
space valued functions (see, e.g., [H]): 

Proposition 2.1. Every u e W^''^{0,T;X) can be redefined on a subset M C 
[0,T] of zero measure so that u G C{0,T]X). (Here and further C{Q^T; X) 
denotes the space of continuous functions /: [0,T] X.) 



Remark 2.2. Let u denote the redefined version of u. Then G C°°([0,T]), 
Va, 6 G [0, T] 



for a.e. t G [0,T]. 

Let X be a reflexive Banach space (with dual space X*) and let iJ be a 
Hilbert space for which there exists a linear bounded dense embedding n: X ^ 
H . Let tt: H ^ X* be the embedding given by tt: h {h, k{-))h, where (•, ■)h 
is the dot product in H. Then embeddings tt and t = t: o k are linear, bounded 
and dense. Triple {X, H, X*) (with embeddings k, tt, i) is said to be an evolution 
triple [11] . For given evolution triple let 



This space is referred to as Sobolev-Lions space T2!|. It is a reflexive Banach 
space with norm given by 



V: -{p,div), 
A: u.^-((u,.)), 
div : u i-> — (u, V-). 





W^^\0,T;X)^{f \ feL^{0,T-X), i{f) e W'^^{0,T; X*)} . 



11/11 



\\f\\L^O,T-X) + \Hf)t\\mo,T;X')- 
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Embedding i is often omitted and the space W^''^{0, T; X) is then introduced as 
the space of functions belonging to L^{0,T; X) whose weak derivative belongs 
to i^(0, T; X*). In this paper i, k and tt will be omitted when they are not the 
subject matter. 

The introduced space has the following property (see, e.g., [11 ): 

Proposition 2.3. Every u € W^''^(Q,T]X) can he redefined on a subset Ad C 
[0, T] of zero measure so that u S C(0, T; H). 

Remark 2.4. Let u denote the redefined version of u. Then Va, 6 e [0, T], Vu G X 
and V(/7 e C°°([0,T]) 

b l-b 

{ut,v) tpdt = {u{t)ip{t),v)H - {u,v)H'Ptdt. 



Consequently, ifyv G X and V<p G C§"{[0,T)) 

{ut,v) ipdt = -{uo'f{0),v)H - / {u,v)H^tdt, 

Jo 

where uq G H, then u{0) — uq. (Similar procedure can be used to identify the 
initial value of a function from VF^'^(0, T; H).) 

Remark 2.5. The mapping 1 1— > ||w(t)|||^ is absolutely continuous with 

dt\\u{t)\\'H = 2{ut{t),u{t)) 

for a.e. i G [0,T]. 

In this paper we consider evolution triples {Hl{Df ,L'^{D)^ ,H-^{Df) {k G 
N) and {V{D), J{D), V{D)*). For both of them the embeddings k, tt and % are 
given by k: u i-> u (natural embedding), tt: u i— ^ /d('^' ') ^^'^ i = tt o k. 

The following theorem describes the space which is dual to a Lebesgue- 
Bochner space [TT] : 

Proposition 2.6. If X is a reflexive Banach space, then 

LP{0,T;X)* = LP'{0,T;X*), 

where 1 < p < oo, l/p+ 1/p' — 1 and duality is given by {f,g) = fgdt, 
f eLP'{0,T;X*), geLPiO,T;X). 

Different constants which are not dependent on the principal parameters 
(such as initial conditions) will be denoted by the same letter C. The dependence 
of such constant on some parameter will be indicated in the subscript. 

2.5 Auxiliary Inequalities 

Let us recall two well-known statements: 
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Proposition 2.7. Let a > 0, 6 > and J be real numbers. If < a + bJ then 
J < 6 + Va. 

Proposition 2.8 (Gronwall's inequality). Let I be an absolutely continuous 
nonnegative function of a variable t £ [0,T] and let Lp,ip L^{0,T) be nonneg- 
ative functions. If the derivative I'{t) satisfies 

I'{t) < (p{t)I{t) +ip{t) a.e. on [0,T], 

then for a.e. t € [0, T] 



I{t) < eJ'o '^(^)''^ (^I{0) + ^(t) dr^ 



The proof of Proposition 12 . 71 is elementary and the proof of Proposition 12.81 
can be found e.g. in [13]. 

We will use the following mix of Lemmas 12.71 and 12.81 

Lemma 2.9. Let I and J be absolutely continuous nonnegative functions of 
a variable t e [0,T], J G L'^{0,T). Let a,c € L^{0,T) and b € L^{0,T) be 
nonnegative functions. If for a.e. tE [0,T] 

I'{t) + J2(i) < a{t)I{t) + b{t)J{t) + c{t) (2.1) 

then 



||J||l2(0,T) < Ca (VW) + ^Ml^O.T) + I|6||l2(0,T)) , 

(/(0) + ||c|| 

where constant Ca depends only on A = ||a||/^i(o.T) 
Proof From for a.e. t £ [0,T] 

I'{t) < a{t)I{t) + b{t)J{t) + c{t) 
and then, by Proposition 12.81 and Cauchy-Bunyakovsky inequality, 



I{t) < exp adr 



/(O) + / {bJ + c)dT ] < 



< eii«"-^(o.*) (/(o) + \\bh2^o,t)\\J\\mo,t) + Mmo.t)) 

and hence 

L2(o,T) + ||c||l1(0,T) + ^(0)) . 

Integrating the inequality < al + bJ + c — I' and noting that I{T) > we 
obtain 

/ J'^dt< [ aIdT+ [ bJdT+ [ cdr + /(0)< 
Jo Jo Jo Jo 

< l|-^l|L~(0,T)||a||Li(0.?') + \\mL-2{Q-T)\\J\\L2{0.T) + ||c||l1(o,T) + ^(0) < 

< Ca (||6|U2(o,t)PIIl2(o,t) + I|c||li(o,t) + 1(0)) (2.2) 
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where 

a = l + ||a|Ui(o,T)e"''"-^'°.-'. 
Applying Proposition 12 . 71 to (|2.2I) we get 

II^IIl^(0,T) < Ca (y^||c|Ui(o,T) + VW) + IIH1l2(0,T)) 

Finally, by Young's inequality 

<elHl.i(o,T, (||6||i.(o,j,) + ||c|Ui(o,T)+/(0) + 
+ ^C^a(li6|liw) + ll^lli^(o.^^)+^(0)) 

= ^a(||fo|li2(o,T) + l|c||Li(0,T)+/(0)), 

where 

=ell»ILi(o,T, j^l + ^C^^ . □ 
2.6 Properties of Special Subspaces 

Let G{D) denote the orthogonal complement of J{D) in L^(D)''. Let Pj and Pq 
denote the orthogonal projectors of L'^{DY onto J{D) and G{D) respectively. 
Projector Pj is referred to as Leray projector. 
Let 

V{D)^ - { /| / e i^^'P)', /(u) = 0, Vu e V{D)} . 

It is clear that functionals represented by Vp, p £ L^{D), belong to V{D)-^. To 
show that any functional from the latter space can be represented is such form, 
the following statements are needed (see fT31 and [H I, §1]): 

Proposition 2.10. A functional f e H^^{D)'^ is reprcscntahlc in the form 
for some p G L'^{D) if and only if 

(/,u)=0 Vuei(7?). 

Proposition 2.11. If p e L'^{D) then 

where the constant C depends only on the domain D . 



|L°°(0,T) 



< e 



l'^lltl(0,T) 
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Propositions l2 ■ lOl and l2 . 1 f] allow us to introduce an operator, which is inverse 
to V: 

Lemma 2.12. The operator V : L'^{D) V{D)^ has a hounded inverse V . 

Proof. For / g V{D)^ let pi — p — Jj^pdx, where p is given by Proposition 
12.101 If there exists p2 € L'^(D) such that / — Vp2 then V(pi — P2) = and 
from Proposition [5TTT] pi — p2 = 0, therefore we can introduce the operator 
V""'^ f ^ Pi (which is clearly linear). 

It follows from Proposition l2.11l that ||pi||£2(£)-) < C'II/II//-i(_d) and therefore 
V^^ is bounded. □ 

Later we will use the operator which was first constructed by M.E. Bogovskii 

Proposition 2.13. Suppose that the domain D is star-shaped with respect to 
some ball. Then there exists a bounded linear operator B: L^{D) — )• Hq{D)'^ 
such that V/ G L^{D) the field u = B{f) satisfies 

div u = / (a.e. in D). 

3 Linearized Equations of Low Compressible Fluid 
Motion 

Consider viscous barotropic fluid with state equation p = Fiji), where p and p 
are density and pressure respectively. Let us linearize this equation near some 
reference pressure Pref ■ 

p ^ F(jp) « F{pref) + F'{pref){p ~ Pref)- 

For brevity assume that p^ef ~ and denote F{0) ~ pq > 0, F'{0) ~ a. From 
the physical point of view l/a=^ = c^>0, where c is the speed of sound. 
Ultimately we have p = po + ap. 

We are going to study the following equations: 

Pt + podivu = 0, 
poUt + Vjj = /lAu + TyVdiv u + pf , (3.1) 
p = Po + ap- 

Here 

p = p{x,t) (p(a;,i)eM), p = p{x,t) {p{x,t)eR), 

xeDci{o,T),te (0,r) c M, 

u ~ u{x, t) is the velocity (u(a;, t) G R''), 

{ — f{x,t) is the external force (per unit volume). 



8 



The coefficients of viscosity ^ > and 77 > are constant throughout the paper. 
The constant a > is referred to as the compressibility factor |^ . 

Equations (13.11) arise as the hnearisation of the equations of compressible 
fluid motion near some solution {d, q} to the corresponding equations of incom- 
pressible fluid. The convective terms {v,'V)ap and (v, V)u are omitted for the 
sake of technical simplicity. 

Let us consider the following initial-boundary conditions for the equations 

m-- 

u\qd = 0, u|t=o = uo, p\t=o = Pq- (3.2) 
It is clear that the solution to the problem p.ip . p.2p depends on the compress- 
ibility factor a: 

Our goal is to study the convergence of these solutions when compressibility 
tends to zero, i.e. a — >• 0. But before doing this, we need to state (and prove) 
the existence and uniqueness theorems for the problem p.l|) . p.2|) and for the 
corresponding incompressible system. 

3.1 Existence and Uniqueness of Weak Solutions 

Consider a non-homogeneous form of the problem (|3.ip . (13. 2^ : 

pt + podiv u = CT, p^ap, 

poUt + Vp = pAu + r^Vdiv u + pf + s, 

u|aD=0, u|i^o = Uo, p\t=o^Po- (3.4) 

Setting cr = and s ~ pgi we clearly obtain the problem p.ip . p. 21) . 
Assume that 

a = L^{0,T;L'^{D)), s e L^{0,T; H'\Df), 

ieL°^iQT)\ n^eL^Df, poeL\D). 

Definition 3.1. A pair 

{u,p} € L\Q,T-HI{dY) X L\0,T;L\D)) 

is called a weak solution to the problem (|3.3I) . p.4p . if for all ip e {[0,T)), 
V G H^{DY and r € L^{D) 

{p,r)^ptdt - {apo,r)ip{0) + / {podiv u, r) ip dt = / {a,r)ipdt, (3.5) 



Po(u, v) (^f - Po(uo, v)(p(0) - / {p,divv)pdt = 

Jo 

= —p / {{u,v)) ipdt — ij / {div u, div -v) If dt + 
Jo Jo 

+ [ {pi,v)^dt+ [ {s,v)ipdt, (3.6) 
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where p — ap. (See subsection 12.21 for the notation.) 
Remark 3.2. From (|3.5I) and (13.61) it foUows that 

peW^'^{0,T;L^{D)) 

and the equations p.3p hold for a.e. t e [0,r] in sense of notation, introduced 
in subsections 12.31 and 12.41 Consequently, for a.e. t G [0,T] 

{pt + podivu,p) ^ {a,p), 
po (u4, u) — {p, div u) = (/iAu + r/Vdiv u + pf + s, u) 

Due to Propositions 12.11 and 12.31 the functions u and p can be redefined so 
that they have well-defined values u{t) and p{t) at each t e [0,r]. Integrating 
the equations above with respect to t using Remarks 12.21 and 12.51 we obtain the 
energy equality: 



Up,\uf + ^\p\'- 
2 V Po 



+ {p\\uf + T]\divu\^)dt = 

= / {api,u)dt+ [ {s,u) dt+ [ —{a,p)dt, (3.7) 
J? Ji Po 

where e [0,T] are arbitrary. 

By Remark |2.4[ from equations (|3.5p . (13. 6p we conclude that u(0) = Uq and 
P(0) =Po- 

Definition 13.11 is equivalent to the following one: 
Definition 3.3. A pair 

{u,p} e L^{0,T;H^{Df) x L\Q,T-L\D)) 



is called a weak solution to the problem (|3.3p . (|3.4I) . if the equations p.3p hold 
in sense of distributions on Qy, and the initial values of u and p are Uq and Pq 
respectively. 

Proof of equivalence of \3.1\ and \S.S\ To prove the implication l3.3l =» IXT] consider 
V e Hl{Df and r £ L^{D). There exist (v„) C C^{Df and (r„) C C^(i:') 
such that v„ — > V and r„ ^ r in Hq{D)'^ and L^{D) respectively. Then, for 
each ^eCg°{0,T) 

*„(x, t) = Vn{x)ip{t) e I?(Qt)^ 
$„(x,0 ='^n(a;)(^(i) e2?(0T). 

Integrating by parts the distributional formulation of the equations p.3p and 
passing to the limit when n — >■ cxd we obtain p.Sp and p.6p with (f e C^(0, T). 
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Then, from Propositions 12.11 and 12.31 it follows that p.Sp and p.6p hold with 
arbitrary e C^([0, T)). 

The implication 13.11 ^ 13.31 follows from the fact that the set 



N 



k=l 



TV e N, (^fe G 2?(0, T), rfe G 2?(D), fc = l..iV 



is dense in 'D{Qt) (see e.g. [S]). □ 

Now we are ready to state the main result of this section: 

Theorem 3.4. The problem (j3.3l) . (j3.4l) /las a unique weak solution {u,p}. For 
this solution the following estimates are valid: 

l|u||L2(o,T:Hi(_D)<i) + l|u||Loo(o,T;L2(£i)<i) + 

+ Va||p||L==(0.T;L2(D)) <Cf^aE, (3.8) 

ll"llw'i.2(0,T;Hi(D)<') - ^'^'^^ 

where 

E = |luo|lL2(£,)d + \fa\\pa\\L'^(D) + 

+ ■y^l|o-|lL2(0,T;L2(Li)) + |ls|lL2(o^T;_f/-i(£')<i) 

and the constants Cf^a,Cf^a depend only on \\{\\l°°(Qt)'' '^^'^ For fixed f 
i/iese constants are bounded when a — )■ 0. Moreover, if po G L^(D) and a G 
L^{0,T;L^{D)), then 

pG M^i'2(o,T;L2(i:))). 

Proof. We will use Faedo-Galerkin method. Let {e^}"^]^ be an orthonormal 
with respect to ((•, •)) basis in Hq{DY and let {ej}°^i be an orthonormal basis 
in L^{D). Then, since Hq{DY is dense in L^{D), for m G N there exist 

m 

i=i 

m 

u° - V c^^^e ■ c^™^ G R 

J = l 

such that p'^ — >■ Po in L'^{D) and — Uq in L'^{DY when m ^ oo. 
Let us look for the solution in the form 



Um = Y^ {t)ej , Pm = X! 4"^ (^)^ 
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where the functions Cj^\t) and c^™\t) are absolutely continuous. Consider the 
following "approximate" system: 

{Pmj. + Podiv u™, a) = (cr, e^) a.e. on [0, T], 
{poUm,t + ^Pm-, e.j) = (mAu™ + ?7Vdiv u™ + p„f + s, e^) a.e. on [0, T], 
i = 1..TO, = apm, 

This system represents a Cauchy problem for a linear system of first order 
differential equations. Under our assumptions for each m it has a unique solution 
{c5™\t),c5'"^(0}™i such that 

Since u„j(t) e span{ej}"^j and Pm{t) G span{ej}JL;^, the following equa- 
tions hold a.e. on [0,T]: 

(Pm,t + /5odivu,„,p„) = (cr,Pm), 
(poUm,t, Urn) - (p™, div u„j) = (pAu„i + TyVdiv u,„ + p^f + s, u„j) 

Hence, recalling the notation from subsection 12.21 by Remarks 12.21 and 12.51 

^ '^Po|u™p + — -|-Ai||u„|p + 77|divu™p = 

= (apmf, u„) + (s,u„j) + — (cr,pm). (3.10) 

Po 



2 V Po 



Let us make some auxiliary estimates. First, from Young's inequality. 



and {a,pm) < ^ (a|p„ip + i|crp). Second, (s,u„i) < ||s||^-i(£,)d ||u,„|| . Then, 
denoting 

I{t) = 1 (^P0\U„M^ + ' = VMllUmWII, 

a = 1 + V^||f||i<»(Q^)d, -^||s(i)||^-i(£,)d, = ^^|(T(t)p 

from p. 101) we obtain 

I' + <al + bj + c 
From this inequality, by Lemma 12.91 we have 

l|J|U^(O.T) <Cf,„(|u^| + V^|p^| + 

+ —^\\'^\\l^{0,T;L^{D)) + I1s||l2(0,T;H-1(£')'')^ ' 
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|/|U~(o,T)<Cf,„(l<J^ + ab°,P + 

1 

+ - 

a 



''■|Il2(o,T;L2(D)) + l|s|li2(0,T;H-i(D)<') 



From the proof of Lemma [^31 it follows that the constant Cf_Q depends only on 
l|f IIl°°(Qt)'' ^^'^ ^1 ^^'^ '^^^ each f this constant is bounded when a — > 0. 

Thus we have proved that estimate (|3.8p holds for the "approximate solu- 
tions" Um, Vm- Since any weak solution to the problem p.3p . p.4p satisfies the 
energy equality p.7p . the arguments above imply that p.Sp holds for each weak 
solution to the problem (13.31) . p.4p as well. 

By Alaoglu-Bourbaki theorem (and Proposition 12. 6p . there exists a subse- 
quence (for brevity not renumbered) such that 

u,„-u inL2(o,r;iji(i5)rf), 
u,„^u inL°°(0,T;L2(i?)'*), 
V^^V mL°"iO,T;L\D)) 

and the limits u and p satisfy (13.81) . 

Now let if e C^([0,r)) and for r e ^2(1?) denote 

Fjn{r) = / (-p„<^t 4- (podivu™ - (t)(^, r) - (ap™(0),r).^(0), 
Jo 

It is clear that (Fm) is bounded in L?{D)* and for each r G spanjejlJL^ F„i{r) = 
0. Then, since {ej}°^]^ is an orthonormal basis of L^{D), 

lim Fra{r) = 0. 

m— )-oo 

On the other hand, 

lim F,„(r) = F(r) ee / {-pipt + (podiv u - a)(p, r) dt - {apo, r)ip{0), 



where p — ap. Thus we have shown that p.Sp holds. Similar arguments for the 
sequence of the functionals 

Gm(v) = / {-po\irn^t,^)dt- I (p„ , divv).^ - po ("m (0) , v).^(0) - 
Jo Jo 

-/ (^Au™ -f TjVdiv -f p^f -I- s, v) di, w e Hl{D f 
Jo 

imply that p.6p holds. 

From p.Sp and p.6p we conclude that a.e. on [0,r] 

\pt\ < -(|podivu|-H|(T|), 
a 

l|utlk-i(r>)<J < II - Vp + /iAu + r/Vu + apf + s\\H-i{Dy / Pa, 
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hence u satisfies (I3.9p . 

If the problem (|3.3p . p.4p had two weak solutions {u^,p^} and {u^,p^}, 
then the difference {u^ — u^,p^ — p^} would be a weak solution to the problem 
p. 31) . p.4p with zero Uq, po> cr and s. Hence, from estimate p.8p we would have 

— u'^ = and — =0. Consequently, the problem p.3p . (|3.4p has only 
one weak solution. 

Finally, if po e L'^{D) and a e L'^{0,T;P{D)), then taking {ej}°^i as an 
orthogonal basis of L'^{D) we obtain p e W^''^{0, T; L?{D)). □ 

Corollary 3.5. Suppose that 



Then the problem p.ip . p.2p /las a unique weak solution {u,p}. 
3.2 Incompressible Limit 

When a = 0, equations p.2p formally turn into nonsteady Stokes equations for 
incompressible fluid: 



Consider the following initial-boundary conditions for the equations p. lip : 



The problem p. lip . (I3.12p is well-studied, so we will only state some stan- 
dard results for it. 

Definition 3.6. A pair 



is called a weak solution to the problem p. lip , p. 121) . if p. lip holds in V^Qt) 
and for aU Lp e C^{[Q,T)) and v e V{D) 



Proposition 3.7. Let f e L'^{0,T; H-^iD)"^), € J{D). Then the problem 
p. 111) . p.l2p has a unique weak solution such that 



feL°°(QT)^ UoeL^iDf, PoeL^iD). 



diYV — 0, 
PoVt -t- V9 = pAv + Pot. 



(3.11) 



(3.12) 



{v,q}eL'{0,T;V{D))xV'{QT) 




V e W^'^{0,T;V{D)), 
q = dtQ, QeC{0,T;L\D)). 
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Moreover, if Vq e V{D) and f G L'^{0,T; L'^{D)), then 

VteL\0,T;J{D)), 

q^L\Q,T-L\D)). 

Remark 3.8. The first part of tlie proposition is proved in [S]. Tlie case Vq S 
V{D) is treated by a minor modification of tlie proof from [F. 

Remark 3.9. Every weak solution to the problem p. lip . (|3.12l) satisfies the 
following energy equality: 

po\v\% + 2iij^ \\vfdt = 2j^ ipo{,v)dt (3.14) 

for ah e [0,T]. 

Regularity of the solutions to the problem (I3.11|) , p.l2p has been rigorously 
studied by V.A. Solonnikov. Let us state a corollary of one of his results (see 
m P- 126]): 

Proposition 3.10 (V.A. Solonnikov). Let £ be a nonnegative integer number. 
Let dD e C2^+2, f e W^f '^(Qt)"* and e J{D) n W^f +2 (£>)'*. Assume that 
all the necessary compatibility conditions up to order i are satisfied. Then 

Remark 3.11. It follows from Proposition [330] that Vq G W^^'^{Q,T; L'^{n)). 
Applying V-i to Vq we get q G W'-^{Q,T;L'^{^)). Then for ^ > 1 

g|t=o- V-iPG(At^o + f(0)), 
i.e. the pressure q has a well-defined initial value. 

4 Convergence Towards the Incompressible Limit 

In this section we suppose that the assumptions of CoroUarv 13.51 are satisfied. 
Let {uq,Pq}o<q<i denote the family of the weak solutions to the problem (|3.ip . 
(|3.2p . Let {v,q] denote a weak solution to the problem p. lip . p.l2p . 

4.1 Convergence of the Velocity 
Theorem 4.1. //a — > 0, then 

u^^v inL'{Q,T-Hl{DY), (4.1) 

xia^v m L°°{0,T;L^{Df), (4.2) 

Vp^^Vg inH-\QTY, (4.3) 
where {v,q} is the solution to (j3.1ip . p.l2p mf/i 

t>o = PjUq. (4.4) 
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Proof. Consider the equality p.6p with arbitrary v G V(rt): 

po{ua,\-)^ptdt - po{uo,v)ip{0) = -p / ((uq, v)) + 

Jo 

+ / (ap„f,v)(^dt + / {po{,v)ipdt, (4.5) 
Jo -/o 

where (/s G C^([0;T)). Clearly (uo,v) = (PjUo,v). The estimate im- 
plies that (uq) is bounded in L'^{0,T; Hq{D)'^) and (yapo!) is bounded in 
L°°{0,T; L^{D)). Then, by Alaoglu-Bourbaki theorem, for any sequence a„ — >■ 
0, n — > oo, there exist w G L'^{0,T; Hq{DY) and a subsequence (nfe)^^^ such 
that 

u,„^ mL^{0,T;H^{Df) 
when fc — > oo. Passing to the limit in (|4.5p . we obtain 

-/ Po{w,v)iptdt- po{PjUo,v)ip{0) = 
Jo 

= -fi[ {{w,v))ipdt + 0+ [ {po{,v)ipdt, (4.6) 

JQ Jo 

which also holds (by continuity) for each (p G Hq{Q,T). Passing to the limit in 
(|3?5|) assuming that r G P{D) and ip G C^(0,r) we get 



/ (div w, r)ipdt = 0. 
Jo 



Since divw G L'^{0,T; L^{D)), this yields that w G L'^{0,T;V{D)). 

Let us collect all the terms of (|4.6p on its left-hand side and denote the 
resulting expression by $(v, (p). For any fixed v G Hq{D)'^ 

$(v,.) Gi/-i(0,T). 
Integrating (|4.6p by parts we see that <i>(v, ip) = 0(v, t)ipt dt, where 

H'^it) = -Poiw,v) ~ p. {{w,v))dt+ / {pof,v)dt 
Jo Jo 

Clearly (/){-, t) G C{0,T; H-\dY). From dHJ and Lemma [2ll2] it follows that 
there exists a unique Q G C(0, T; L'^{D)) such that for i G [0, T] 

Let p denote the distributional derivative dtQ in V^Qt)- Then is 
a weak solution to the problem (|3.1ip . (|3.12l) with the initial condition (|4.4p . 
By Proposition 13.71 this solution is unique, so {w,p} — {v, q}. Hence the whole 
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family converges to v when a — ^ 0, i.e. (|4.ip holds. (Otherwise this family 
would have a subsequence which wouldn't converge to v weakly. But such a 
subsequence would be bounded and applying Alaoglu-Bourbaki theorem and 
the arguments above we would come to a conclusion that Uq converges to v, 
which is a contradiction.) 

As a consequence of estimate (|3.8p and Alaoglu-Bourbaki theorem, we can 
pass, if necessary, to another subsequence of (rife), so that u^^^ ^ w in L°°(0, T; L'^{DY) 
for some w. Hence V^s € C^(0, T) and Vr G L'^iDf 

I {ua,^^,r)Lpdt ^ / {•w,r)ipdt, fc — > oo. 

From this equation and (|4.ip we conclude that w — v for a.c. t e [0,T]. Thus 
Ua„^ V in L°°{0, T; L^{D)'^). Since the whole family converges to v when 
a -> 0, reasoning as above we conclude that (|4.2[) holds. 

Finally, (14. 3p follows from (|4.ip and the passage to the limit in p.6p when 
a^O. □ 

Theorem 14. II is in good agreement with the results proved in [5J[T]. But the 
topology of the convergence (|4.ip , (|4.3I) can be strengthened using the approach 
coming from the artificial compressibility method JF: 

Theorem 4.2. Suppose that Vq = PjUq. When a — > 0, 

Ua^v in L^{0,T;H^{Df), (4.7) 

Ua-^v in L°°{0,T;L^{Df), (4.8) 

Vp^^Vq inH-\QTY, (4.9) 

if and only if Uq G J{D). 

Proof. By Remark El and Theorem O Vv e V{D) 



Po(Ua(t) - Uo, v) + 



-/I / ((Uq, v))(it - 7/ / (div UQ)(divv) (it + / {pa + ap){^,\r)dt 
Jo Jq 

^ -fJ- iiv,v))dt-0+ po{f,v)dt^pQ{v{t)-Vo,v). 



By estimate (|3.8p |uQ(t)| is bounded for < a < 1. Then, since V{D) is dense 
in J{D), by Banach-Steinhaus theorem Vv e J{D) 

(u„(t),v) ^ (tJ(t),v). (4.10) 

Now consider a family of real numbers 

x„ = po|u„(r) - tj(r)|2 + -|pjT)|2 + 

Po 

+ 2^ / \\ua - vW^ dt + 2tj / |div(u„ - r')|^dt. 
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Expanding the parentheses we write 
where 



=_2po(u„(T),tJ(T))-4/x / {{uc„v))dt, 

Jo 

l-T i-T 

X'i^ = po\n^{T)\^ + -\p^{T)\^ + 2n \\uJ^dt + 2Tj \divu^f dt. 
Po Jo Jo 

By energy equahty p.l4l) 

= Pol-^ol' + / {poi,v)dt. 
Jo 

From (PIUI) and dUI]) 

^ -2p,{v{T)MT))-^P f{{v,v))dt = -2X(i). 

Jo 

Now we pass to the hmit in energy equahty p.7p using estimate p.Sp : 
Po Jo 

^Po|uoP+0+ / (pof,'y)di = ^*'^+Po(|uop-|t;on. 

"'0 



Hence we have proved that po(|uoP — ji'oP) when a -> 0. Then (|4.7p and 

(|4.8p hold if and only if |uo| = l^ol = |PjUo|, i.e. Uq G J{D)- And the passage 
to the hmit in ([5^ implies (g^. □ 

4.2 Convergence of the Pressure 

The important difference between the equations of compressible fluid and the 
equations of incompressible fluid is that only the former enjoy the mass conser- 
vation property. Namely, let 



D 



M{t) = I pdx^ I {po + ap)dx. 

Then ([XT]) implies 



dM 
-;- = 0. 

dt 

Consequently ^j^p{t)dx doesn't depend on t. In contrast, this is not true for 
the incompressible fluid, since q can be shifted by arbitrary function of t. 
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However, for arbitrary A £ K the pressure q can be shifted so that q{t) dx — 
A for a.e. t £ [0,T]. (Let C{t) = (/ q{t)dx - A)/ dx, then the substitution 
q{t) := q{t) — C{t) yields the desired result.) 

Theorem 4.3. Suppose that Vq ^ Uo E J{D) md, 

/ q{t)dx^ / podx for a.e. t e [0,T] (4.11) 
Jd Jd 

Then 

Pa ^ q in L'^{0,T;L^{D)), a 0. 

Proof. Assumptions of the theorem imply that the solution {v, q} to the "in- 
compressible" problem (|3.1ip . p.l2p belongs to the same functional space as 
the solution to the "compressible" problem (|3.ip . (|3.2p . Then the difference of 
the solutions {Uq,Pq,} = {uq, — v,pa — g} is a weak solution to the problem 
I, p.4p with initial conditions 



U„|t=o = 0, P^lt^o = Po = po - q{0) 
and with non-homogeneous terms given by 

a = —aqt, s = aqt. 
Assumption implies qt G L'^{0,T;L^{D)). Then ([331) implies 

\\'^a\\L^{Q^T-Hi{D)'') + \\^a\\L°-={0,T;L^{Dy) < Cy/aE, 
II^"IIl°°(0,T;Z2(d)) < CE, 
E = ||-Po!|£2(£,) + ||9t||L2(o^T;L2(D)) + Va||f||oo||'7||L2(0,T;L2(n)) 

thus we have proved the convergence of the velocity (once again). 

Uniform boundedness of the pressure Pq, by Alaoglu-Bourbaki theorem, 
implies that for every sequence a„ — >■ 0, n — t- oo, there exists a subsequence 
{nk)ken such that 

Pa„, in L°°(0,T;L2(£)))^ fc oo, 

where P* G L°°{Q,T;L'^{D)). Passing to the limit in the equation (|3.6p (as in 
the proof of Theorem 14. ip we show that VP* is zero. Then P* = and the 
whole family P^ converges to zero. □ 

Remark 4.4. The estimate p.Sp not only yields the fact of convergence, but 
also gives the rate of the convergence. In other words, (13. 8p allows to estimate 
the error of approximation of a compressible fluid by the incompressible one. 

It would be desirable to have strong convergence of the pressure. But in 
general case such convergence is impossible, since (due to Proposition [23J 

- '7||l~(0,T;L2(i5)) > IIpo - g(0)||L2(D)- 
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Theorem 4.5. Suppose the assumptions of Theorem \4--^ are satisfied and the 
domain D is star-shaped with respect to some ball. If, in addition, 

qeW^-^{0,T;L^{D)) and Po^q{0), 

then 

Pa-^q in L°°{0,T;L^{D)), a ^ 0. 
Proof. Let w = —B{qt)/pQ. By Proposition l2.13l 

and 

l|w||n/i,2(o,T;//i(D)<i) < C\\qt\\w^-^(0,T:L^{D))- 

Then {Uq, Pq} = {u^ — v — aw^pa — g} is a weak solution to the problem p.3p . 
p.4p with initial conditions 

Vc.\t=o - Uo = -aw{0), P„|t=o = Po = Po - 9(0) = 

and with non-homogeneous terms given by 

cr = 0, s = a(— poWt + yuAw + ?7Vdiv w) + aqf. 

Clearly 

\\^\\L-2{0,T;H-i-{D)<i) < aCCi, 

where 

Cl = l|w||i4^l,2(o,r;ffi(£l)<i) + l|f l|oo||9||L2(0,r;L2(£))). 

Assumption ((4ll|) implies qt E L'^{0,T;L'^{D)). Then, from jSH) 

l|Ua||i2(o,T;ifi(-D)'') - 

ll^a|lL°°(0,T:Z2(D)) < C'-E^, 

where 

-E^ = (l|Uo||L2(n)d + V"ll^b||L2(i5) + 

which yields the desired convergence when a — > 0. □ 
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